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1 INTRODUCTION 



ABSTRACT 

We present a critical analysis of the usual interpretation of the multicolour disc model 
parameters for black hole candidates in terms of the inner radius and temperature 
of the accretion disc. Using a self-consistent model for the radiative transfer and the 
vertical temperature structure in a Shakura-Sunyaev disc, we simulate the observed 
disc spectra, taking into account doppler blurring and gravitational redshift, and fit 
them with multicolour models. We show not only that such a model systematically 
underestimates the value of the inner disc radius, but that when the accretion rate 
and/or the energy dissipated in the corona are allowed to change the inner edge of the 
disc, as inferred from the multicolour model, appears to move even when it is in fact 
fixed at the innermost stable orbit. 

Key words: accretion, accretion discs - black hole physics - radiative transfer - 
X-rays: stars 



sion from the accretion disc, not only is the disc assumed to 



X-ray *pectra ot galactic black liole candidates (UtitlUj 



be geomet rically thin and optically thick ( shakura & Sun 



are becoming an increasingly important tool in determin- 
ing both the physical properties of the source surroundings 
and, as the ultimate goal, the parameters (mass and spe- 
cific angular momentum) of the black hole. Two main spec- 
tral states, as defined by their spectral components and flux 
level (typically in the 1 — 10 keV band), have been observed 
in GBHC. In the Hard/ Low state, the sources emit most 
of their power in a hard power law tail with photon index 
r ~ 1.3 — 1.7 and exponential cutoff at about 100 keV. In 
the generally more luminous Soft/High state, most of the 
energy comes from a thermal component with characteris- 
tic temperature kT < 1 keV, while a power law compo- 
nent (with r ~ 2.0 — 2.5) dominates above a few keV. In 
addition, sources have been reported to be in a Very high 
state, spectrally intermediate between the Soft/High and the 
Hard/Low states, when the power law component has a flux 
comparable to the thermal one and no sign of a cutoff is seen 
in the high energy tail. 

We will focus our discussion on the Soft/High state. 
According to the common paradigm in interpreting the 
observed spectra, the ultra-soft thermal component is the 
product of the emission from an optically thick accretion 
disk. On the other hand, the hard power law component is 
probably due to multiple Compton scattering of soft photons 
by a population of hot electrons which reside in an active 
(or flaring) coronal region surrounding the accretion disc. 

In the zeroth order approximation for describing emis- 



yaev 1973), but also the vertical temperature structure is 



neglected, as is the effect of Comptonization on the emer- 
gent spectrum. Each point of the disc is then assumed to 
radiate like a blackbody at an effective temperature which 
scales with the radius as r~^^^ . This is the so-called multi- 
colour disc model (MOD; Mitsuda et al. 1984). Such a simple 
model has the advantage of being easy to use in trying to 
fit spectral data. It has only two adjustable parameters: a 
[colour) temperature, independent on the source distance, 
and a normalization factor, which in turn depends on the 
inner radius of the disc, on the distance of the source and 
on the inclination of the disc from the line of sight. When 
these two latter quantities are independently known, or can 
be reasonably inferred, fitting the ultra-soft component of 
a GBHC spectrum with the MCD model immediately gives 
estimates of the temperature in the inner region of the disc 
and its inner radius. 

The aim of this letter is to demonstrate that the val- 
ues one can obtain with this fitting procedure, which is now 
a standard one, are not in general directly related (or at 
least not in a reliable way) to the actual disc parameters. 
We will use a self-consistent model for the radiative trans- 
fer in Shakura-Sunyaev accretion discs about a compact 
source, including first-order corrections for gravitational red- 
shift and Doppler blurring, as developed by Ross & Fabian 
( 1996| ). The observable spectrum is computed for a num- 
ber of different physical situations, each determined by the 
values of the viscosity parameter a, the accretion rate and 
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t he fraction of the accreted power dissipated in the corona 
(Svensson & Zdziarski 1994). In all the cases we consider, 



the value of the inner radius of the disc will be kept fixed. 
Nevertheless, when trying to fit the computed spectra with 
the MCD model, we will obtain results at variance with this 
assumption. That is because changes in the normalization 
factor of the MCD model are produced, in a very compli- 
cated way, by the variations in the disc structure (and hence 
in the Comptonized spectrum) induced by variability of the 
accretion rate and coronal activity. 



2 STRUCTURE OF THE DISC 

We assume the basic structure for spatially thin accretion 
discs around Schwarzschild black hole s given by the standard 
theory of Shakura & Sunyaev ( 1973| ). The disc is assumed 
to have a fixed inner boundary at the innermost stable orbit 
at radius iiin = 3i?s — &GM/c?. We consider both an in- 
ner, radiation pressure dominated region and an outer, gas 
pressure dominated region of the disc. In determining the 
structure of the disc, the opacity is assumed to be domi- 
nated by Thomson scattering. 



0.2(1 +X) cm^g-\ 



(1) 



where X is the mass fr action of hydrogen. Following Svens- 
son & Zdziarski (1994), we slightly modify the standard set 
of disc equations, allowing a fraction / of the disc accretion 
power to be dissipated in the corona instead of in the cold 
disc itself. Furthermore, at every given radius R, the den- 
sity is taken to be uniform in the vertical direction (which 
is the correct solution for the disc structure when radiation 
pressure dominates). 

Choosing dimensionless parameters 
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where MEdd ~ 3.1 x 10~*m Mq yr~^ is the Eddington ac- 
cretion rate for a disk efficiency rj = 0.083, the flux emerging 
from the surface of the disc is given by 



Fo = 1.2 X 10^''m"V~^[mJ(r-)](l - /), 



(2) 



with J(r) = (1 — •y/3/r). The radial structure of the disc is 
summarized by the values of the half-thickness of the disc 
(in uni ts of the Schwarzschild radius) h and the uniform gas 
density po: 



(i) Radiation pressure dominated region 
h= ^ = 5.3(1+ X)[mJ{r)]{l~f) 

o A . -I n-7 -1 -1 3/2 

Po = 3.4 X 10 a m r ' 
x[mJ(r)]-^(l-/)-3gcm-^ 

(ii) Gas pressure dominated region 

/I = 3.4 X 10-2(1 -f X)^/^0a-^/^V-^''^°r2^''2° 

x[mJ(r)]^/^(l-/)^/"' 

po = 4(1 + x)-^/^«a-^/i«m-^/i°r-3^/20 

x[mJ{r)f{l - fr'^^^^g cm-^ 
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The radiative transfer and the vertical temperature 
structure are treated self-consi stentl y at each radius, as 
described in Ross & Fabian (1996), using the Fokkor- 



Planck/diffusion equation of Ross, Weaver & McCray (197J ) 
in plane-parallel geometry. The local temperature profile, 
T{z), is found by balancing the heating rate due to dynamic 
heating, Compton scattering and free-free absorption with 
the cooling rate due to inverse Compton scattering and free- 
free emission. 

Incoherent Compton scattering within the accretion 
disc must be treated properly whenever it has an important 
effect on the emergent spectrum. The competition between 
Compton scattering of higher-energy photons and inverse 
Compton scattering of lower energy photons often results in 
the emergence of a Wien-law hard tail. The effective optical 
depth for absorption is given by 



(7) 



where tt is the Thomson depth below the surface, and Tg{i') 
is the optical depth due to free-free absorption. For high en- 
ergy photons the thermalization depth (where Ttiy) = 1) can 
be reached at very high Thomson depths, where the temper- 
ature is considerably higher than near the surface. Compton 
scattering in the outer layers of the disc downscatters these 
photons to lower energies and produces a Wien-law tail. 

Decreasing the flux emerging from the disc itself, either 
by lowering the accretion rate (m) or increasing the fraction 
(/) of the accretion power dissipated in the corona, results 
in ever larger regions of very high density, both in the outer- 
most portions of the disc and, for very low fiuxes, in the very 
innermost portions as well. In those cases Comptonization 
is not complete or 'saturated' because 



TTt < 1 



(8) 



at the thermalization depth. In these regions we treat Comp- 
ton scattering as coherent, dropping the Fokker-Planck term 
in the radiative transfer equation (see discussion in Ross & 
Fabian ( |l99e| )). 

In our calculations we fix the value of the black hole 
mass m = 10 an d assu me a composition X — 0.71. Following 
Ross & Fabian (1996) we divide the region 3 < r < 200 of 



the disc into 20 annuli that make comparable contributions 
to the total luminosity. For each set of the parameters a, m 
and /, we find the radius at which radiation pres sure equals 
gas pressure, w hich is the root of the equation ( ^hakura & 



Sunyaev 1973) 



TAB 



o^nl^ \2/21- 16/21 /, J•^6/7 

370(ma) ' m ' (1 — /) '. 



(9) 



J(rAs)i«/2i 

We use equations (3) and (4) for all the annuli for which 
3 < r < rAB, and equations (5) and (6) for r > rAB- The 
typical emergent spectrum from each annulus is calculated 
and then multiplied by the area of the annulus to find the 
contribution to the spectral luminosity. The resulting spec- 
tra are added together, taking into account blurring due 
to gravitational redshift and transverse Doppler effect us- 
i ng th e method described by Chen, Halpern & Filippenko 
( [1989D . Finally, in order to allow direct comparison with ob- 



servations in which a dramatic change in the inner radius 
has bee n reported (GRS 1915-1-105, see e.g. Belloni et al. 
(1997t)), we choose an inclination angle i = 70° and a dis- 
tance D — 12.5kpc. 
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3 THE MULTICOLOUR DISC MODEL 

For each set of parameters the emergent disc spectrum is 
used as a table model to create a fake set of data (with 
RXTE response matrix p012_LRl_970804.rsp, integration 
time 1000 s). These data are then fitted with the standard 
muhicolour blackbody model (DISKBB in XSPEC) in the 
2 — 20 keV band. The model assumes that the local emis- 
sion from the disc is Planckian, with a temperature profile 
T(r) cx r~^^^. Therefore the observed flux from the disc is 
approximated by 



Table 1. Physical parameters. 
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where B(E, T) is the Planck function and Tout is the tem- 
perature at the outer radius of the disc (and it is assumed 
that Tout ^ Tcoi). The two fit parameters of the model are 
the colour temperature of the inner accretion disc (Tcoi) in 
keV and the normalization factor 



R. 



]J2 ' 

where _Rcoi is expressed in km and the distance D in units 
of lOkpc. 

The parameter Rcoi is not the effective inner disc radius 
(i.e. the radius at which the temperature of the disc is the 
highest). The assumption behind using the MCD model is 
that the Comptonized emergent spectrum can be approxi- 
mated by a diluted blackbody spectrum with a colour tem- 
perature (Tcoi) higher than the effective temperature (Toir). 
The ratio fmi = Tcni/Tcff , the spectral hardening factor 



( Shimura fc Takahara 1995 ), is then assumed to be constant 
throughout the disc and for varying physical parameters (rh, 
a, f, etc.). Thus the actual inner edge of the disc would be 
given by 



Rin = vRcs = m{i)Rcoi{Tcoi/T,«Y 



(11) 



where r) is the ratio of the inner radius of the disc to 
the radius at which the emissivity actual l y peaks, typically 
77 ~ 0.6 - 0.7 ( IShakura fc Sunyaev 1973| ; [shimura fc Taka-' 



hara 19951 ). The factor [T^oi/T^sf comes from the assumed 



dilution of the blackbody spectrum, while g{i) takes into 
account gener al relativistic correc tions and is of the order 
3 ~ 0.7 - 0.8 ( lEbisawa et al. 1994| ). 

In the next section we will show that these assumptions 
are not justifled in general due to the differences in Compton 
processes induced by changes in m, /, and a. 



4 RESULTS 

We have simulated disc spectra in eleven different physical 
situations represented by the ten sets of parameters Sl-Sll 
(see Table |^) . The sets are listed in order of decreasing value 
of the radius tab, the boundary between the radiation pres- 
sure and gas pressure dominated regions. Table |l| also lists 
the value of the radial coordinate rcoh for the boundary of 
the region(s) where Compton scattering has been treated as 
coherent (see section ^). As the power dissipated in the disc 
decreases and the density increases, from the combined ef- 
fects of changes in m and /, this boundary moves inwards. 
That is, for an increasing portion of the disc Comptoniza- 
tion is incomplete and the local spectrum will be somewhat 
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0.5 
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0.1 
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13 
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0.1 


0.05 


0.8 







Table 2. Fit parameters in the range 2-20 keV. All the models 
refer to a IOMq black hole with inclination angle and distance set 
to the values of the black hole candidate GRS 1915-1-105, namely 



i = 70° 


and D = 


12.5kpc. 
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56.9 


22 


1.76 


S3 


0.75 


53.3 


6.0 


1.82 


S4 


0.75 
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S6 


0.69 
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46.5 


2.8 


1.95 


S8 


0.68 


41.1 


1.6 


2.08 


S9 


0.67 


41.5 


1.5 


2.07 


SIO 


0.59 


39.0 


1.8 


2.12 


Sll 


0.55 


24.5 


0.7 


2.68 



harder. For the SIO model, the density in the innermost re- 
gion of the disc is so high (due to the J(r-)~^ term) that we 
have to consider incomplete Compton scattering also in the 
first inner annuli, where a significant fraction of the X-ray 
power is emitted. Consequently, for this set of parameters, 
two values of rcoh are listed. Finally, for the model Sll, 
when the disc is entirely gas pressure dominated, we consid- 
ered incomplete Compton scattering throughout all the disc 
(but see next section for a caveat on the applicability of our 
model to this extreme case). 

Figure 1 shows the observed spectra for SI, S4, S7 and 
SIO, taking into account Doppler blurring and gravitational 
redshift. For comparison, Fig. 1 also shows the blurred spec- 
tra that would be observed if each point on the surface of 
the accretion disc emitted a blackbody spectrum at the lo- 
cal effective temperature. The observed spectra are harder 
than such effective blackbody spectra, and the discrepancy 
increases as we increase the density of the disc and decrease 
the power released in the disc itself (going from SI to SIO). 

Table ^ lists the values of the fit parameters Tcoi and 
J?coi for the ten cases, as well as the corresponding reduced 
of the fits and the spectral hardening factors obtained 
from eq. ( pT| ) by assuming rjg{70°) = 0.5 and Rin = 88.6km, 
as appropriate for the IOMq Schwarzschild black hole we 
are considering. The results are also shown in Figures 2a 
and 2b, where we plot the spectral hardening factor and the 
colour radius as functions of m(l — /), which is a measure 
of the flux emergent from the disc itself. 

To check the sensitivity of our results to the energy 
range fitted, we also fitted the spectra in the 3 — 10 keV 
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Figure 1. The disc spectra (solid lines) computed for the parameter sets SI {m = 0.3, f = 0, a = 0.1); S4 (m = 0.2, / = 0.5, a = 0.3); 
S7 (m = 0.1, / = 0.4, a = 0.1) and SIO (m = 0.05, / = 0.5, a = 0.1). For each of the plots also shown are the spectra that would be 
observed for a blackbody-emitting disc (dashed lines). 



band. The values obtained for the colour radii were, at most, 
just one or two kilometers smaller than those obtained in the 
2-20 keV band (i.e. a change of about 3 — 4%). 

It is clear that the multicolour model gives systemati- 
cally lower values for the "inner radius" than the actual ones, 
a point which is well known. More interestingly, our results 
show that the hardening factor is by no means constant when 
we change the physical parameters of the disc. For example, 
assumi ng the conventional value /coi = 1.7 ± 0.2 (Shimura 
&i Taka hara 199f^ ; Sobczak et al. 199£| ) and taking Doppler 
blurring and gravitational redshift into account, we would 
have underestimated the value for the inner disc radius by 
up to a factor of 2 (Sll) using the fits with a multicolour 
model. On the other hand the model gives quite stable and 
acceptable results for high accretion rates and/or lower val- 
ues of the fraction of the power eventually dissipated in the 
corona. This in turn can help to understand the cases in 
which a nea rly constant value of _R coi is observed even when 
Tcoi varies (Tanaka & Lewin 1995). 



5 DISCUSSION 

We have shown that MCD models systematically underesti- 
mate the value of the inner radius of the accretion disc for a 
black hole candidate. We also have shown that the spectral 
hardening factor /coi, which is needed to correct the results 
of the fits with a multicolour model, is not, as is usually as- 
sumed, constant when the accretion rate and/or fractional 
coronal activity change. 

Recent observations of galactic black hole candidates 
seem to point towards an extreme dependence of the ob- 
served colour radius on /. 

In the work of Sobczak et al. (1999), where the au- 



thors report on the RXTE spectral observations of the 1996- 
97 outbursts of the microquasar GRO J1655-40, a sudden 
jump inward of the color radius (see their Fig. 7) occur for 
the five very high state observations of the source during 
which the power-law (hard) fiux was extremely high and 
the blackbody-to-total fiux ratio was less than 0.5 (see their 
Fig. 5). 



A similar trend has been reported by Muno et al. ( 199t ) 



in the case of the microquasar GRS 1915-1-105. In their ob- 
servations, the smallest values of the inner disc radius ob- 
tained by fitting the spectrum with the MCD model are 
associated with the highest values of the power-law-to- 
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cal structure of the disc, its density profile and, in particular, 
the surface density. 

It is quite possible, and perhaps easier to conceive, that 
dramatic changes in the disc spectrum are produced by 
changes in the surface properties of the accretion disc rather 
than by the disappearance of the entire inner regions. This, 
unfortunately, is hard to model theoretically. It would be in- 
teresting, for example, to assess the problem in the extreme 
case when almost all the accretion power is dissipated in the 
corona (i.e. in the limit / ^ 1). As we believe that even in 
our most extreme simulation (Sll) with / = 0.8 the simple 
treatment of the vertical disc structure we used here can be 
inadequate, and as this must be true a fortiori for larger 
values of /, we do not push our study further. 

To do that we would need to take into account the X-ray 
reflection spectra of discs illuminated by a hot corona and 
to model carefully t he hea ting of the disc atmos phere (e.g., 
see Sincell & KroUk ( |l997| ), Ross & Fabian ( |l993| ), who show 
how a hard tail in the disc spectrum is produced in this case) . 
It is natural to believe, given the trend we observe in the 
simulations presented above and the observational results 
we referred to, that in such a situation the oversimplified 
multicolour model will give even smaller values of the colour 
radius. 

A detailed treatment of this kind requires further nu- 
merical work which is beyond the aim of this Letter. 



Figure 2. (a): the points are the values or Hcol (in km) plotted 
as functions of m(l — /), which is a measure of the flux emergent 
from the disc, while the solid line marks the actual inner boundary 
of the disc, fixed in all our ten simulations, (b): triangles represent 
the value of the spectral hardening factor needed to rescale the 
computed iJcol to Ri^. 
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blackbody flux ratio, which in turn should be directly related 
to the value of /. 

In both these cases the inferred inner disc radii can 
change by more than a factor of four. 

Belloni et al. (1997a,1997b), analyzing RXTE data of 
GRS 1915-1-105, interpreted the sudden change measured in 
the inner colour radius of roughly a factor of four in terms 
of the rapid disappearence of the whole inner part of the 
geometrically thin accretion disc. Once again, in that case 
the inner edge of the disc appears to shrink to its smallest 
values during an outburst in which the power law flux dom- 
inates the blackbody, and the blackbody-to-total flux ratio 
is less than 0.5. 

This behaviour is exactly what our results predict we 
should expect from the MCD model when we increase / 
and/or reduce the accretion rate. In this case, every time 
we the observations imply that the coronal activity is domi- 
nant, the multicolour fits should be corrected with a varying 
hardening factor 1.7 < /coi < 3 in order to recover the actual 
value of the inner disc radius. 

It should be clear, however, that the issue of determin- 
ing the exact shape of accretion disc spectra is very com- 
plicated and is very sensitive to the actual physical models 
used to describe it. In particular, key elements are the verti- 
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